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Abstract 

We formulate an adiabatic theorem adapted to models that present an 
instantaneous eigenvalue experiencing an infinite number of crossings with the 
rest of the spectrum. We give an upper bound on the leading correction terms 
with respect to the adiabatic limit. The result requires only differentiability 
of the considered spectral projector, and some geometric hypothesis on the 
local behaviour of the eigenvalues at the crossings. 
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1 Introduction 



The availability of intense pulsed laser sources has opened a large field of possibilities 
to control atomic and molecular dynamical processes. One of the main theoretical 
tools to analyze these processes is adiabatic Floquet theory [|GMDJ|| and references 



therein. The general setup can be described as follows. One considers a molecule 
described by a Hamiltonian Hq acting on a Hilbert space 7i, in interaction with 
one radiation mode of frequency u. (The description of the interaction with several 
modes of different frequencies can be formulated along similar lines). Since the 
intensity of the field is quite large, the field is treated as a classical field. The 
Hamiltonian of the molecule perturbed by the electromagnetic field can be written 
for example as 

H = H + EM F(ujt + e ) (1) 

where M is the dipole moment of the molecule and E e K is the amplitude of the 
radiation field, F is a 27r-periodic function and 9q the initial phase. We assume that 
Hq has a discrete spectrum. In order to describe a laser pulse the amplitude is taken 
as a slowly varying time dependent function E(et), where one takes e.g. e = 1/T p 
with T p the duration of the pulse. A new technique that provides an efficient method 
for complete transfer of population is based on frequency chirping: within the pulse 
duration the frequency is also slowly modulated lu = uj{et). 

This model has thus two kinds of time-dependencies in the Hamiltonian: one 
that is periodic and another one that is slowly varying. The periodic part can 
be treated by Floquet methods, and the slowly varying part by adiabatic theory. 
Adiabatic Floquet theory is based on the following statement: Assume that in the 
Hamiltonian ([I]) the parameter E and the frequency u are made time dependent, 
E(t), u>(t). Consider the propagator U (t, to; 9 ), solution of the Schrodinger equation 

% ^U(t, t ; 60) = H{ut + 6 ) U(t, t ; 9 ), U(t, t; 0„) = 1 (2) 

acting on the Hilbert space 7i. We consider an enlarged Hilbert space by tensoring Ti 
with the space of square integrable functions on the unit circle: L 2 (§ 1 ,7i). The op- 
erator U(t, to; 0) can be lifted into the enlarged space, interpreting the ^-dependence 
as a multiplication operator. We can then define 

U K (t, t ) = e" M * )e U(t, t ; 6) e toU> ^ d 

where d = 3/ 36. The statement is that equation (pD is equivalent to 

i^U K (t,t )=K(t)U K (t,t Q ) (3) 

with 

K(t) = - tw (t)-^ + H + E(t) M F(6) 

and zu(t) denotes an effective instantaneous frequency defined by w(t) = u(t) + 
tdui{t)/dt. Assuming that the time dependence of E(t), u(t) is slow one can develop 
adiabatic techniques for the evolution of When K has pure point spectrum, the 
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first ingredient are the instantaneous eigenvalues and eigenvectors. They always can 
be written and labelled in the form 



Xj^k = Xj ; o + tu k, k <E Z (4) 

The index j has the same cardinality as the dimension of the Hilbert space TC. Thus, 
even if we take simple models with finite dimensional 7i, the Floquet spectrum has 
infinitely many eigenvalues. As functions of E and w, these eigenvalues may exhibit 
crossings, which the adiabatic approximation can accomodate in case there is a 
finite number of them, see |BF ], | |Hag| | . The structure (Q) of the eigenvalues is 



such that if we consider a slowly varying frequency w(t) that goes through 0, one 
can encounter situations in which a branch of instantaneous eigenvalues undergoes 
an infinite number of crossings with other branches. Let us stress here that this 
situation is not generic, as actual crossings of eigenvalues are more the exception 
than the rule. However, we give below a whole class of systems for which this 
situation is true. Note also that in case w(t) passes through 0, the domain of 
K(t) becomes time dependent, so that technical issues regarding regularity of the 
evolution operator have to be addressed. This is done in the Appendix. 

The goal of the present paper is to formulate an adiabatic theorem that can 
be applied to such situations with an estimate on the corrections to the adiabatic 
limit. Adiabatic Theorems without gap conditions are known to be true, see [AE|, 
however, in general, no estimates on the error terms are available. 

While this work was motivated by the physical situation described above and 
discussed below in the examples, our analysis of the adiabatic approximation is 
model independent and can be applied to more general situations. 



2 Adiabatic Theorem 



2.1 Context 

The adiabatic approximation in Quantum Mechanics has a long history which we 
will not attempt to retrace here. We refer the reader to the recent surveys |AE2|, 



| JPfl and references therein. Let us simply recall here that the works following that of 
Born and Fock |BF| by Kato ||Kal|| , Nenciu ||Nen|| and Avron, Seiler, Yaffe |[ASY|| have 
lead to a formulation of the Adiabatic Theorem under the usual gap assuption that 
is general and where the error term is well controlled and of order e. In case the gap 
assumption is modified, the situation is less explicit. In this section, we switch back 
to the notation H(et) for the slowly varying time-dependent Hamiltonian. Assume 
H(s) is smooth in s G [0, 1] and there exists a spectral projector P(s) of H (s) which 
is strongly C 2 on [0,1]. Avron and Elgart have shown in [|AE|| that the adiabatic 
theorem holds under these conditions, provided P{s) is of finite rank, independently 
of any spectral considerations. The limitation of this approach is that, in general, 
no estimate can be made on the rate at which the adiabatic regime is attained. In 
certain specific situations, an estimate on this rate is available. In the case where 
the spectral measure fi v is a-H61der continuous, with = P'(s)ip, ip being the initial 
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condition, the rate of convergence was shown in ||AE|| to be of order e a /( 2+a ) . A case 



where the spectrum of H(s) is assumed to be dense pure point is dealt with in [ |AHS| . 



Another situation, considered in [|Joyj| , where the gap hypothesis is not necessarily 



fulfilled occurs when H(s) = H (s) + eHi(s), where the domain of Hi(s) is smaller 
than that of H (s). In both cases, the error term remains of order e. In the present 
article, we consider another situation in which the usual gap assumption is modified 
and the error made in the adiabatic approximation can be estimated. We make 
the hypothesis that the spectral projector P(s) is associated with an eigenvalue 
A(s), in the sense that H(s)P(s) = X(s)P(s), for all s G [0,1]. We assume that 
A(s) is isolated in the spectrum except at a series of times {ofej^gN accumulating 
at a g]0, 1[ where it experiences crossings with the rest of the spectrum. Requiring 
some conditions on the local behaviour of the gap between A(s) and the rest of the 
spectrum near the crossing points Ok, we estimate the error term in the theorem 
without a priori knowledge on the nature of the rest of the spectrum. Note that for 
s = Ok such that A(ofc) is not isolated in the spectrum, P(ok) does not represent the 
entire spectral projector associated with the eigenvalue A(ofc). 

2.2 One crossing 

Let us make more precise the regularity hypotheses under which we shall work. In 
order to deal with the application described above, we will assume the hamiltonian 
is unbounded. This causes technical difficulties motivating the part ii) of the hy- 
pothesis below which justifies our manipulations. We show in the appendix that 
this assumption is verified for our models. In case H(s) is bounded, this part of the 
assumption is automatically verified. 

HO) i) We assume that for all s G [0, l]\{a}, H(s) is a strongly C 1 self-adjoint 
operator defined on a dense domain T> independent of s in a separable Hilbert 
space /C, where < a < 1. Whereas if (a) is bounded self-adjoint on /C. We 
also assume the existence of a spectral projector P(s) of H(s) which is strongly 
C 2 on [0, 1] and such that H(s)P(s) = P(s)H(s) = X(s)P(s), for all s G [0, 1]. 

ii) Further assume that the unitary evolution operators U(s) = U (s, 0) and 
A(s) = A(s,0) generated by H(s), respectively H(s) + ei[P'(s), P(s)} (see (|), 
(|D) are well defined for all s G [0, 1] and possess the properties i) to v) listed 
in Theorem |A7 



Note that P(s) needs not be finite dimensional and A is continuous. 
We start by considering one crossing of A with the rest of the spectrum by 
revisiting the strategy proposed in [IBFJ , making use of the general analysis presented 



in 



Let g(s) be the gap between A(s) and the rest of the spectrum of H(s): g(s) = 
dist (A(s), <r(s)\{A(s)}) > 0, s G [0,1]. We also introduce the bounded, strongly 
C 1 operator L(s) = i [P'(s), P(s)]. We assume that g^jO} = {o} and consider the 
strong differential equations on T> 

ieU'(s) = H(s)U(s) , U(0) = 1 (5) 
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ieA'(s) = (H(s)+eL(s))A(s), A(0) = 1. (6) 

The unitary A is the so called adiabatic evolution which possesses the well known 
intertwining relation A(s)P(0) = P(s)A(s) ||Kal|| , |[Kr|| . Finally, let W(s) be defined 
by W(s)= A- 1 (s)U(s). We have on V 

iW'(s) = -A- 1 (s)L(s)A(s) W(s), W(0) = 1, (7) 

in the strong sense. To compare the adiabatic and actual evolutions, we need to 
compute the size of the difference of the unitary W(s) at two times surrounding the 
crossing. This is the aim of the next result. 

Lemma 2.1 Under the above assumptions, we have for any 0<uo<t<o<s< 
ui < 1, 

|| W(uq) - W{u x ) || < C (e\u - t\/g 2 t + e\u x - s\/g 2 s + e/g t + s/g s + \s-t\) (8) 
where g t = inf ue[uo . 

t] 9{ u ) > 9s — i n fu£[s,«i] 9{ u ) an d the constant C is uniform in uq, 

u\, s and t. 




Figure 1: The various quantities defined in Lemma pZTT 



Remark: 

On the basis of the classical paper by Born and Fock, 



|BF|| , and the detailed analysis 



of crossings by Hagedorn (Hag], one would expect the corresponding estimate with- 
out the first two terms. However, such an estimate requires more detailed knowledge 
of the structure of spectrum, e.g. that the gap is given by the distance between two 
eigenvalues, than what we assume in our general setting. 



Lemma 2.1 can be used to treat two standard situations: 



1) If there is a gap G between A(s) and the rest of the spectrum, this lemma 
implies that the adiabatic approximation holds with and error term bounded 
by Ce/G 2 . 
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2) If one starts the evolution on a crossing point which splits like s a near 0, we 
can use this lemma to show that the adiabatic approximation is valid with an 
error bounded by 

\\U{1)- A{1)\\ <Ce 1/{1+2a) 

if e is small enough. This is precisely the situation encountered at the begin- 
ning of the interaction of a laser pulse with frequency that is in resonance with 
the difference between two energy levels of the molecule [ |Ho]| , |GJ|| . 

To get this estimate, we can consider only half of the problem by letting aside 
all the terms containing a t and setting u\ — 1: 

\\W{1) - W{Q)\\ <C (e\l-s\/g 2 s + e/g s + s) (9) 

This is indeed fully justified by the proof of the lemma (see below). Next, 
we have by hypothesis that g(s) > g s = G s a if s is small. Introducing this 
behaviour in Equation (|), we obtain \\W(1) - W(0)\\ < C (e/s 2a + s). The 
result follows now by balancing the two contributions by choosing s = s(e) = 
£ i/(i+2a) _ Ag a fn ; with more information on the spectrum, as in ||BF| , [Hag 
should be able to improve the above estimate to order eVCH-a). 



one 



Proof of Lemma |2.1| : 

The idea of the proof is to integrate Equation (|7|) over the interval [uq, ui] and then 
to get "nice" estimates of the sizes on each subintervals [tto,£], [t,s] and [s, By 
integrating Equation (^), we get 

i {Wiux) - W(u )) = - f A- l (u)L{u)A{u) W{u) du (10) 
A-\u)L{u)A{u) W(u) du - r A- 1 (n)L(?i)A(u) W{u) du 



For the middle term, we simply use the properties of the operator norm and the fact 
that A(u) and W(u) are unitary to obtain: 

\\W(s) - W(t)\\ < f \\L(u)\\ du < sup \\L(u)\\ \s-t\ , (11) 

h «e [o,i] 

i.e. we do not care about the behaviour of g(u) inside the subinterval [t, s]. 

To estimate the first integral, let Q(u) = 11 — P(u). A simple computation, using 

P(s) P'(s) P(s) = 0, shows that 

P(u) L{u) P{u) = Q{u) L{u) Q{u) = , (12) 

and due to the intertwining property of Aiu), we can write 

W(t) - Wiu ) = i f (p(0)A- 1 (m)L(m)A(m)Q(0) 

Juo ^ 

+ Q(0)A- 1 (w)L(n)A(M)P(0))iy(M) du (13) 
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Now, we need to extract an explicit e dependence from this equality in order to 
obtain the estimates stated in the lemma. To do this, we follow [ |AS Y|| and introduce 
the bounded operator 1Zl{u) defined by 



TlAu) = — i R(u, X)L(u)R(u, A) dX 
2m Jr(u) 



» 

where R(u, A) = (tl(u) — A) is the resolvent of H (u) at A and where the loop 
r(w) is a circle centered at A(u) of radius g{u)/2. It has the properties (see [|ASY|| , 

[K L {u),H(u)) = [L(u),P(u)} (14) 
P{u)K L {u)P{u) = Q{u)K L {u)Q{u) = 0. (15) 

Standard arguments show that TZl{u) is strongly C l and that 

KM = 7T- i f i?(it, X)L'(u)R(u, A) 
2m Jr(u) K 

- R(u, X)H'(u)R(u, X)L(u)R(u, A) (16) 

- R(u, X)L(u)R(u, X)H'(u)R(u, X))dX , 

where H'(u)R(u, A) is to be understood as the bounded operator 

H'(u)R(u, A) = H'(u)R(u, i) (1 + (A - i)R(u, A)). (17) 
Hence, we get the following estimates: 

\\n L {u)\\ < ^^\\L{u)\\{g{u)/2)- 2 = 2\\L{u)\\/g{u). (18) 

271 

\\K(u)\\ < cmax{\\H'(u)R(u,i)\\\\L(u)l\\L'(u)\\}/g 2 (u). (19) 



The main property of TZl(u) (see |[ASY|| ) is that it satisfies for any ip e V the 



following equalities, as verified by means of (|i~4"|): 

P(Q)A-Hu)L(u)A(u)Qti))il) = -ie 4- (P(0)A-\u)K L (u)A(u)Q(0)Tp 

nil V 



du 

+ ie P{G)A-^{u)K' L {u)A{u)Q(G) x/j . (20) 



and 



Q(0)A- l (u)L(u)A(u)P(0) ^ = is— (Q(Q)A-\u)TL L (u)A{u)P{<d) if; 

- ie Q(0)A- 1 (m)^' l (m)A(m)P(0) . (21) 

These equations imply that /* A^ 1 (u) L{u)A{u)W (u)du is proportional to e. Indeed, 
Equalities (|1^) and the intertwining property of A(u) show that the diagonal blocks 
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are 0. 

Introducing Equalities and (PTj) in Equation (|T3"D, we get 



W(t) - W(«o) = (22) 

— £ 



t ^(Q(0)A- 1 (w)^ i (w)A(w)P(0)-P(0)A- 1 (M)^(«)A(M)g(0))w(w)dw 



j du 

-e ( (P{0)A-\u)n' L {u)A{u)Q{0)-Q{0)A-\u)n' L {u)A{u)P{0))w{u)d 



Performing an integration by part in the first integral, using the differential equation 
(0) for W(u) and taking into account that A(u), W(u) are unitary and P(0), Q(0) are 
projectors, gives us the following bound for the norm of the difference W(t) — W(uq): 

\\W{t) - W(uo)\\ < 2e (\\n L (t)\\ + \\TZ L (u )\\ 

+ sup \\K L (u)\\\\L(u)\\(t-u ) + sup \\K' L {u)\\{t-u )). (23) 

ue[uo,t] «e[«o,t] 



Next, we use first Estimates (£jj|) and (|T^) and then the fact that < uq < t < 1 to 
obtain the desired bound: 

\\W(t)-W(u )\\ < - sup ||L(«)|| + - sup \\L(u)\\ 2 (t-u ) 

9t «e[«oi*] 9t ue[u ,t] 

+c^| sup {\\H\u)R(u,i)\\\\L(u)l\\LXu)\\\(t-u Q ) 

9t we [wo,*] 

< 12- sup {||L(u)||,||L(u)|| 2 } (24) 

#t «e[o,i] L J 

+2 £^^ {11^)^^)1111^)11^1^)11} 

9t «G[0,1] 

/£ e\t-u \ 

< c 2 — + 



0t & 2 



Using the same kind of arguments, shows that on the subinterval [s,iii], we have 



\\W( Ul ) - W(s)\\ < c 2 (tl^li + 1\ . (25) 

9s 9s 

Combining Estimates fllTD, (Q) and ( 23 ) gives the announced bound for ||14 7 (m 1 ) — 
W(u Q )\\. □ 



2.3 Infinite number of crossings 

We now have all the information required to proceed to the case of an infinite 
number of crossings. We make the following hypotheses describing what happens in 
the neighbourhood of each crossing. 

Spectral Hypotheses: There exist two partitions {w^j-fcgN of [0, a) and (a, 1] re- 
spectively: 

= < • • • < u k-l < U k ■ ■ ■ * U ^o = a = u£ < • • • «fc < Ufc-l <■■■<< = 1 

such that for each A; e N*, 



S 



HI) one can find non empty open intervals V k , which satisfy V k C [u k _ 1 ,u k ], 
K C [ut, ut- x ] and 

sup flf(s) < inf </(*) (26) 

where 4" = u k }\V k and # = [u^ , u^]^ . 

H2) there are constants G±(k) > and a fc-independent positive constant a such 
that for all s G V k : 

G ± (k)\s-o±\ a <g(s), (27) 

for some points o k G . 




Figure 2: Illustration of the spectral hypotheses H1-H2 on the interval (0, a). The 
intervals Vp are represented by ( ). 

Comments: 

-/ ^ These Spectral Hypotheses mean that the crossings are well separated and that 
they behave as power of order at most a. Hypothesis HI tells us that outside the 
crossing regions (V k ) the gaps are relatively "large". This means that the only 
accumulation point of small gaps is a. 

2) The choice of a constant exponent a is not as restrictive as it might look at 
first. Indeed, we are interested in an upper bound, so it is the greatest a that will 
determine the global behaviour. 

3) In the applications, we will consider examples where g -1 {0} = {o k }\ the set of 
crossing points of A(s) with the rest of the spectrum. This implies a > 0. But, the 
case of an infinite number of avoided crossings can be treated by taking a = in 
Hypothesis H2. 

To obtain an estimate for the difference between the real evolution U(l) and the 
adiabatic one A(l), the idea is to apply Lemma |2.1| on a finite number of cross- 
ings and to take a simple integral bound (as in (JTTJ) ) over the rest of the interval 
surrounding a. The choice of the number of crossings will be optimized with re- 
spect to e in order to get a simple form for the bound of the remainder term. To 
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state the corresponding result, we need to introduce some notations. Let A±(k) = 
max{|u£-o£|, |«fc_i-Ofc |} and r±(Jfe) = max [A ± (k)/G%(k), A%(k)/G±(k)}. The 
functions K \— > \u% — a\/ J2k=i T±(k) l l <yl+2a ^ are monotonically decreasing to zero, 
so, if e is small enough, we define K±(e) G N* as the greatest integer satisfying 



u K — a\ 



1/(1+2q) 



> £ 1 /(1+2 Q ) _ m 



This integer always exists if e is sufficiently small and, by construction, K±(e) — > oo 
as e —* 0. 

Theorem 2.1 For e small enough, under HO and the Spectral Hypotheses HI, H2 
and provided that 

s(sT ± (k)) 1/(1+2a) < \Vf\/2 for all 1 < k < K ± (e) , (29) 

for some constant <; > 0, we have that 

U(l) = A(l) + O (max{|^_ (£) - a\, \u% +(e) - a\}) . 

Hence, as \im e ^ K ± (e) = oo, \\U(1) — ^4(1)|| goes to zero for e — > as fast as 
max{|w^_ (£) -a|,|w+ +(e) -a\}. 

Remarks: 

1) The theorem states that the error can be estimated provided we can compute 
the critical value K±(e). Further considerations on the practical aspects of this 
computation are given in the next section. 

2) Condition ( p9|) implies that the size of the intervals V^r cannot be too small with 
respect to e r±(k). 

3) While we shall apply the theorem in a situation where the spectrum is simple 
and pure point, the theorem remains valid under the sole existence of an eigenvalue 
separated from the rest of the spectrum by gaps with the properties stated in H1-H2, 
without any knowledge on the rest of the spectrum or restriction on the dimension 
of P(s). 

Proof: 

In the sequel, we will denote by the same symbol c all inessential constants. Let us 
consider the interval [0; a). In order to simplify the notations, we will not write the 
sub/super-scripts — . Picking some t, s 6 Vjt such that t < Ok < s and \t — Ok\ = 
| s - o k | , we get 

\\W{u k ) - W(u fc _i)|| < c (e\t - Uk-i\/gt + e\s - u k \/g 2 s + e/g t + e/g s + \t - s\) 
<c{er{k)\t-Ok\- 2a + \t~o k \) (31) 
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by the preceding section. Indeed, we have that g t = inf u6 [ UA ,_ 1)t ] g(u) = g{r t ) for 
some r t G [uk-i,t]. Now, by Hypothesis HI, r t G Whence, we have that 

gt = g{n) > G(k) \r t - o k \ a > G(k) \t - o k \ a 

as rj < t < Of.. Using the same kind of arguments, we can show that g s = 
mf ue [ S}Uh ] g{u) > G{k) \s — o k \ a - Finally to obtain the bound (f30j), it remains to no- 
tice that \s — t\ = \t — o k \ + \ s — o k \ = 2 \t — Ofe| together with \t — Ok\, \t — Uk-i\ < A(fe) 
and \s — Ofc|, \s — < A(fe). 

We now get an estimate by choosing t = t(e, k) in order to balance the two 
contributions appearing in the last term of Equation (|30|) above: for some constant 
<T > 0, we set 

-\t(e,k)-o k \, (32) 



|t(e, k) - o k \ 



2<> 



|t(e,fc)-o fc | =<;(er(k)) 1/{1+2a) . (33) 

By definition, k) G V^, hence, as k will eventually be bounded from above by 
K(e), this imposes Condition (^0) in the statement of the theorem. Replacing t by 
t(e, fc) in (^) and summing over fe, we get for any K < K(e), 

\\W(0)-W(u K )\\<c(, + ,~ 2a )j:(er(k)) 1/{1+2a) . (34) 

k=l 

On the other hand, using the differential Equation (0), we obtain, 

\\W{u K ) -W(a)\\ < [ a \\L(u)\\du<c\u K -a\. (35) 

J UK 

Again, we balance the two right hand sides in ( |34]) and ( j35|) by setting the integer 
K = K(e), which has been defined in Equation (^8|) . Consequently, 

\\W(0)-W(a)\\ < cfc + ^^EH^ + kw-fl 

^ fc=i 

< c (<; + <; _2q + 1) \u K ( e ) -a\= C(q) \u K ^ - a\ (36) 

where C(g) is independent of e. Proceeding similarly on (a, 1] completes the proof. 
□ 

Remarks: 

1) The introduction of an adjustable constant q is necessary in the following appli- 
cation to satisfy the hypothesis of the theorem. 

2) In the step fl3Ti| ) we deliberately lost a little in the estimate by using \t — o k \~ a < 
A£ |t — o k \~ 2a in order to simplify the subsequent arguments. It is nevertheless pos- 
sible to get slightly sharper results by not adopting this simplification, however the 
analysis gets more involved and less transparent. We simply note here that in the 
examples discussed below, this more careful analysis yields, for the generic situation, 
an error term of order e p with an exponent p = 1/3, instead of the value p — 1/3 — u, 
for any v > obtained there. 
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3 Application 



We can obtain more explicit estimates on the rest by considering some specific 
behaviour at the crossings. 

Let us introduce the following notation: F k ~ f(k) means that there exist two 
constants < C\ < c 2 < oo such that C\f(k) < F k < C2f(k) for k G N* large 
enough. We have the 



Proposition 3.1 Assume the hypothesis of Theorem \2. i| and the following behaviour 
for the relevant quantities: 

| w ± - a \ = d/k? + C 2 /k? +1 + o(l/^ +1 )) , (3 > 0, Ci ^ 

G±(k) ~ k\ 

\V±\~l/k 5 , 5>0. 

We set fi = mm{(3+ 1 + 2 7 ,a(/3 + l) + 7 }. T/jen ||C/(1) - A(l)|| = 0{e p ) where the 
exponent p is given by 



P 



1 



1 + 2a 
1 

1 + 2a ' 



if n > (1 + 2a) 
v Mv > if // = (! + 2a) 



< (l + 2a) 



[ (/3 + l)(l + 2a)-/i 
provided that 5 satisfy the following constraints: (3+1 < 5 < /3 + max{l, /z/(l + 2a)}. 
Proof: 

The idea of the proof is to explicit conditions on the different exponents ensuring 
the validity of Theorem [2.1| . We will only consider the interval [0, a), the same kind 
of arguments will apply on (a, 1]. Again, in order to simplify the notations we will 
let aside the sub/super-scripts — . 

First, we have that 2 A(Jfe) = u k -u k ^ = C^/k^ 1 + o(l/k 0+1 ) ~ l/k p+x , which 
implies that 

5 > (3 + 1 > 0, (37) 

since 2 A(k) > \V k \ ~ l/k s . Notice that the length of the V k can be rescaled by a 
uniform constant if 8 = (3 + 1. 

Next, A(k)/G 2 (k) ~ l/k p+1+ ^ and A a (k)/G(k) ~ 1/A^ +1 ) + t. So, if we denote by 
// = min{/? + 1 + 2 7 , a(/3 + 1) + 7} then r(k) = m&x{A(k)/G 2 (k), A a (k)/G(k)} ~ 
1/k^ by increasing the overall constant in Theorem |2.1| if necessary. Whence, 



f>(A:) 1/(1+2a) ~f>-^ 1+2Q ) ~{ log if if/i = l + 2a (3* 



fe=i fe=i 



K° if // > 1 + 2a 

logif if/i = l + 2a 

K l- M /(l+2a) if ^ < 1 + 2« 
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and considering the definition of K(e) (see Eq. we obtain 

' K(e)~P if //> l + 2a 

K(e)- /3 /logK(e) if /i = 1 + 2a 
#( e )-/J-UV(i+2a) if//<l + 2a 



.1/(1+2 



a ) |^( e ) ~ a l 



K{e) 

Y, r(k) 1/{1+2a) 

k=l 



(39) 



Condition ( |29|) stated in Theorem [O] reads 

1/(1+2q) 



< \V k \/2 



(40) 



for all 1 < k < K(e). Notice that this condition is automatically satisfied if 5 < 
/V(l + 2a). In general, it will be satisfied for a sufficiently small q, if 



F(e) = £ l/(l+2«) K ( £ )<5-M/(l+2a) 

remains bounded as e — > 0. Using (|3^), we have 

K ( £ y-f>-r/(i+2a) if /i > l + 2a , 
i J1 (e)~'j ^(e)*- 1_ ^/logif(e) if/i = l + 2a, 
K{e) 5 ^- 1 if/i<l + 2a. 



(41) 



(42) 



As K(e) — > oo for e — > 0, Equation ( f4"2"D implies that ^(e) will remain bounded if 
5 < /3 + max{l;/i/(l + 2a)}. 

Hence, using (|5B|) and (|3"9"D, we get that the Adiabatic Theorem |27I] holds with a 
remainder term on [0, a), 



O (J«A'_( e ) - a l 

where the exponent p is given by 

1 



O K(e) 



V — \ 



1 + 2a 
1 

1 + 2a 



0(e % 



if // > (1 + 2a) 



i/ W > if /i = (1 + 2a) 



[ (J3 + 1) (1 + 2a) - n 



if n < (1 + 2a) 



provided that /3 + 1 < 5 < /3+max{l; /// (l + 2a)}. To determine p in case /i = l + 2a 
and S = (3 + 1, we have used the estimate e-V(i+2«) ^ f^)/ 3 log if(e) < Kie)^' 
for all v ' > 0. This ends the proof of the proposition. □ 



Remark: 

We can strengthen the last comment of the preceding section with the following 
observation. In case a — (3 — 7 = 1 and 5 = 2, this last comment asserts that 
we have p = 1/3, instead of p = 1/3 — u, for all v > 0. Now, if in Lemma 2.1, 
the right member were missing the terms e\u —t\jg\ + e\u\ — s\/g 2 , as one would 
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expect with a little more information on the spectrum, an analysis similar to the 
one provided above leads to an error term of order e 1 ^ 3 . This makes it reasonable to 
expect that in such a situation the error actually is of that order, as it was the case 
in the corresponding analysis of one crossing performed in [BF|, see 
it is shown in the examples below that the values a = (3 = 7 = 1 
generic in some sense. 



Hagfl . Finally, 
and 5 = 2 are 



4 Examples 

We now consider a family of models for which the situation just described takes place 
as the effective frequency w takes the value zero. We start by considering the most 
general model for a two level system driven by a periodic field. The model can be 
characterized by choosing freely the eigenvalues A +im = \ + +mu and A-^ = X- + ku 
and the corresponding eigenf unctions of the form: 

, , ns fe ix ^ cosz(9)\ im0 , . //n f-e~ iy ^ sin z(9)\ ike 

(43) 

in which the functions x, y and z are periodic modulo an integer multiple of 9. 
Defining the unitary matrix 



T(0) 



e ix ^ cosz{9) - e - iy W sinz(6) 
e*»W mn*(0) e~ ix W cos z(9) 



the corresponding Floquet Hamiltonian can be written as (dropping the 9 depen- 
dence in the notation) 

K = -iwd-iwT {dT- 1 ) + T D T L 

where D = diag(A + , A_). Using the notation 2p = x + y, 2$ = y — x and choosing, 
without loss of generality, A + = — A_ = A, the Floquet Hamiltonian can be expressed 
as 

K = -i zud (44) 



+ 



w&d + (A — zudp) cos(2z) ^ — iwdz + (A — wdp) sin(2^)^e 
( %wdz + (A — wdp) sin(2z)^e 2 ^ —wdfl — (A — wdp) cos(2z) 



where df denotes the derivative with respect to 9. 

We will consider two different models with the same eigenvalues but with different 
eigenfunctions. We remark that since the adiabatic theorem depends only on the 
properties of the eigenvalues (and regularity properties of the projectors), it gives 
the same upper bound for the correction for all the models ( fill) with equal spectrum. 
But the actual deviation can be very different from this estimate, depending on the 
efficiency of the couplings. 

We can choose for example the following eigenvalues: 



\± t k(w) = k w ± (r](w) + w)/2 , where r](w) = J (w — lu ) 2 + Q, 2 (45) 
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and uq, Q are constants. The first model is defined by choosing x(9) = 0, y(9) = 9, 
i.e. 2i?(0) = 2p(9) = 9 and cos(2 z) = —{w — ujq)/t]{w), sin(2 z) = fl/r](zu), hence z 
is independent of 9. The corresponding Floquet Hamiltonian is given by 

d , 1 / ujn Vte~ 



which is the usual RWA model widely used in quantum optics. 

The second model is defined by the choice x{9) = —g(9)/2, y{9) = 9 — g(9)/2, i.e. 

2-#(6>) = 9, 2p(9) = 9 — q{9) and the same z as for the RWA case. This leads to 

«=w + ^(r; "ill) 

This Hamiltonian can thus be understood as a perturbation of the RWA model. 

We consider now an extra smooth slow time dependence in the parameter w = 
zu(s) and in g = g(9,s). This implies that the eigenvalues, the eigenvectors and 
the corresponding eigenprojectors are smooth functions of s, so that the regularity 
Hypothesis HO i) is satisfied. We show in appendix that HO ii) is satisfied as well 
for any choice of smooth functions x, y, z and A. 

We assume, for simplicity, that w{s) = s (but any other smooth monotonic func- 
tion of s would equally do). We select the eigenvalue A(s) = A +i0 (s) = (r)(s) + s)/2 
and denote by if) the associated eigenvector. The only crossings that A experiences 
are with the A-^+i's and they take place at times s such that 

r](s)=ks, keZ*. (46) 

We remark however that these crossings can lead to corrections to adiabaticity, or 
not, depending on whether the corresponding eigenvectors are coupled. The non 
adiabatic coupling among the branches is measured by the following scalar product: 

I d s ^ k+1 (s)} = -i- jT 2 Vfc+i)0-2i P (M(y( s ) -i sm(2z( S ))tf'(M) d9 



As) 
2tt Jo 



^ e ike+i e (e,s) de 



where the ' denotes the derivative with respect to s. 

Recall that the couplings between the eigenstate if)(s) associated with the level 
A(s) and its orthogonal complement in the Hilbert space is given by the operator 
L(s) = i[P'(s), P(s)], see (0), since the adiabatic evolution A(s) follows the in- 
stantaneous eigenspaces. A direct computation of the matrix elements (^_jt + i(s) | 
L(s)V ; (s)) with P'(s) =| if)'(s)){if)(s) | + | ^(s))('0 / (s) | shows that the above 
scalar product is proportionnal to the couplings responsible for the non adiabatic 
transitions. 

For the RWA model, as g = the non adiabatic couplings are given by: 

(ip(s) | d s V>-,fc+i( s )} = -z'{s)dk,o 

Thus, the level A(s) is not coupled to the infinitely many other levels it crosses. 
Hence we are lead in this case to an effective problem displaying no crossing, so that 
the error is of order e in this case. 
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Figure 3: The first eigenvalues of the RWA and modified RWA models. 



For the other model, we will obtain non zero couplings at all the crossings, 
if we choose g(6,s) such that exp(i g(9, s)) has infinitely many non zero Fourier 
components. For example, one can take g(9, s) = p(s) sin(#) (in particular p can be 
chosen constant). This coupling is then given by 

(^(s) | d s ^ k+1 (s)) = (-l) k+1 z'( S )J k (p(s)) 

where Jz, is a Bessel function. 



We will now verify that the assumptions of Proposition |3J] are satisfied. Let us 
focus on the interval (0, S], for S small enough. The interval [— S, 0) can be treated 
similarly. Again to simplify the notations we will not explicit the + sub/super- 
scripts. 

In fact the following analysis is valid for all the models fl44|) under the assump- 
tion that the eigenvalues can be written as A± jm (s) = ms ± N(s)/2, where K is a 
C 2 function with bounded derivatives such that K(0) > 0. In particular they are 
satisfied for the eigenvalues given in fl45f ). The Hypotheses imply that the function 
/f (s) = N(s) — C s is strictly decreasing for any ( greater than, say, some Co- Under 
these conditions the following assertion shows that the crossings that A(s) = K(s)/2 
experiences with the rest of the spectrum take place at times such that K(s) = ks, 
k G N large enough. 

Assertion 1 For ( > (o, the function f((s) = — ( s has a unique positive zero 
and if ( < £ we have > o^. 

From the expansion 

/ c ( S ) = K(0) + (K'(0)-C)* + O( S 2 ), 



16 



we obtain the behaviour of o^. 

We define the sequence Uk > by the equation: 

tt(u k )-ku k = (k + l)u k -K(u k ) i.e. N(it fc ) = (k + l/2)u k . (48) 

Assertion [I] implies that u k < Ok < u k -i and, from Equation (refzopos) and the fact 
that Uk = o k+ i/2, we obtain 

u ^ k+ iT*(o) +o{i/k3) - <49) 

Next, we have 

Assertion 2 On the interval [u k , ttfe-i], i/ie spectral gap is given by 

g{s) = dist (A(s), a(s)\{A(s)}) = |N(s) - jfe s| < u fe _i/2 . 

More precisely, for u k < s < Ok we have that g(s) = N(s) — ks < Uk/2 and for 
Ok < s < Mfc-i we /iai>e i/ia£ o(s) = ks — N(s) < u k -i/2. 

This assertion is easily proven by considering the different cases. 
We now prove that the Spectral Hypothesis H1-H2 are verified. Assertion |] and 
Equation ( pEUD show that the sequence is (for k large enough) monotonically 

decreasing to a = 0. To define the intervals V^, we choose any point r k in (o k ,u k -i) 
such that g(r fc ) = kr k — ^(r k ) < u k /2 and set V k = (u k ,r k ). The V k s are disjoint 
and Ik = {uk}U[rk, u k -i). By definition of V k , we have that g(s) < Uk/2 = g(u k ) and 
for r k < s < Uk-\ the gap is given by g(s) = k s — K(s) > Uk/2. Whence, Hypothesis 
HI is satisfied. Finally to prove that H2 holds, we need to estimate the behaviour 
of g(s) on Vk~. the Mean Value Theorem implies that for each s G Vk\{ok}, there is 
an q s , in the interval joining s and Ok, such that 

g(s) = |K(s) - ks\ = \k - #'(q s )\ \s - o k \ ~ k \s - o k \ , 

which shows that H2 is satisfied with a = 1 and G(k) ~ k. 

It remains to check the conditions given in the statement of Proposition |3.1| . We 
have 

\u k - 0| = u k = N(0)/ife + K(0)(K'(0) - l/2)/k 2 + 0(l/k 3 ) i.e. (3 = 1 

G{k) ~ k i.e. 7 = 1, 

[Vfc| ~ 1/k 2 i.e. 5 = 2. 

(50) 

To get the estimate for \Vk\, we have used that (u k , o k ] C V k C (wfe, Wfc-i] and the 
expressions for o k , and in Equations ( $7\ ) and (|45|). This implies that, \i = 
a {j3 + 1) +7 = l + 2a; and 5 = /3 + 1. So, we can use the second case of Proposition 
|3.1| to prove that the adiabatic approximation holds for the models: 

\\U{1) - A{1)\\ < ce p , for any p <- . (51) 

3 

In keeping with the last remark of Section 2, we recall that a more careful analysis 
yields p = 1/3. 
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A Appendix 



In this appendix, we show that an operator on L 2 (5' 1 ,7i) of the form 

K(s,e) = -iw(s)^+H(s,6), (52) 

where H(s, 9) is a bounded operator in TC such that s \— > H(s, 9) and s h-» d/d9H(s, 9) 
are norm continuous and s \— > ccj(s) is continuous, admits a strongly continuous uni- 
tary propagator U(s) = U (s, 0) with all expected regularity properties, even if there 
is a value a for which w{a) = 0. Notice that the assumptions on H will be satisfied 
if, for example, (s, 9) \— > H(s,9) is strongly C 1 . 

The proof relies on a theorem of T. Kato [ Ka2|| , which we will restate in a more 
suitable form for our purpose. 

Theorem A.l (Kato) Let JC and T> be Hilbert spaces such that T> is densely and 
continuously embedded in JC and let K(t), < t < T , be a family of self-adjoint 
operators in /C. Suppose that 

1) T> C domi^) for all < t < T , whence the K(t) are bounded operators from 
D to K,, and the application t \— > K{t) is norm continuous from T> to JC; 

2) there exists a family of isomorphisms S(t) from T> to K, which is strongly 
continuously differentiable and such that 

S{t) K{t) S{ty l = Kit) + B[t) 

where B(t) is a strongly continuous bounded operator on JC. 

Under those conditions, there exists a unique family of unitary operators U{t, s) on 
JC defined for < s, t < T with the following properties: 

i) U(t,s) is strongly continuous on JC in s,t with U(s,s) = 1; 

ii) U(t,r) = U(t,s)U(s,r); 

Hi) U(t,s)V C V, \\U{t,s)\\x> < iVe c l*~ s l and is strongly continuous on V in s,t 
simultaneously; 

iv) —U(t, s)ijj = i U(t, s) K(s) ip for any ^eD, for < s, t < T; 

v) for each ip G V and fixed s, ~^U(t, s) exists and is equal to —i K(t) U(t, s) ip 
and strongly continuous in JC in t. 

To prove this theorem, we apply Theorem 6.1 in | Ka2 | to the operator A(t) = 
iK{t), which is stable with constants of stability c = and N — 1 (see Definition 
3.1 and Theorem 4.1 therein). The fact that U{t,s) is unitary follows from the 
self-adjointness of K(t), the construction of U(t,s) by unitary approximants given 
in the proofs of Theorem 4.1 and 6.1 in |[Ka2|| and the invertibility of U(t, s), which 
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is a consequence of the fact that A°(t) = —iK(T — t) satifies also the hypothesis of 
Theorem 6.1 in ||Ka2|| . See also Remark 5.3 therein. 



We now prove that the family of self-adjoint operators defined by Equation (p2|) 
satisfies the hypothesis of Theorem |A.1| . To simplify the notation, we will not explicit 
the ^-dependence and write d for d/dO. 
Proof: 

For V, we choose dom(— i w*d) for some > 0, and we notice that for any t such 
that vj(t) 7^ 0, we have that dom.fr (t) = D and if w{t) = 0, then dom.fr (t) = /C. 
For the norm on T>, we choose the graph norm associated to —iw*d: 



v=W\\ 2 + \\-iv>*dip\\ 2 > 



Whence, V is a dense continuously embedded subspace of K. 
For any s, t and any ip G T>, we have 

\\{K{t)-K{s))if>\\ 2 < 2 |ro(t) ~f (g)|2 || -iw.di/>\\ 2 + 2 \\(H(t) - H{s)) 




< 2 max 1 " ;\\H{t)-H{s)\\ 2 \ ll^ll^ 2 



which shows the norm continuity of K{t). 

We set S{t) = S = —i w^d + i. S is an isomorphism between V and K which 
is strongly differentiable (by t independence). It remains to show that S satisfies 
Hypothesis 2) of Theorem [A.l| . For this, we first notice that for any if) G dom K(t), 
we have that S 1-1 ^ GPc dom K(t) and 

K{t)S- x il) = S' x K{t)ip + H{t) S^if) - S- x H{t)^ 

= S^Kit) if; + S- X S H{t) S- V - S^ff (t) ip 

= S- 1 (K{t) -iw*dH(t)S- 1 )^ . (53) 

Whence, for any if) G dom K(t), we have that the left-hand side of Equation (|53|) 
belongs to P. So we can write, 

S K(t) S~ V = V - i w*dH{t) S~ l il) , for all ^ G dom £T(t) . 

Setting = —iw*dH{t) S~ x , we have a strongly continuous bounded operator 
(by the assumptions on if) which satisfies SKfyS^ 1 D K{t) + B{t). To show 
the reverse inclusion, we can consider any b > 2 sup t \\B(t)\\ which implies that ib 
belongs to the resolvent set of both K(t) + B(t) and S K(t) S* -1 . It follows that 
{Kit) + B{t) + i b)~ x C S {Kit) + i b)~ S~ x . But since the left hand side has domain 
/C, we must have equality between Kit) + B{t) and S K{t) S -1 instead of inclusion. 
□ 

In the examples of Section 4, both H{s, 9) defined through ( [44]) by means of 
smooth fonctions x, y, z, A of (s, 9), and H {s, 9) + ei[P'{s, 9), P{s, 9)} where P{s, 9) = 
\ift{s, 9)){ip{s, 9) \ with if){s, 9) given by one of the vectors ( fl3"D satisfy the hypotheses 
of the Theorem. Hence assumption HO ii) is satisfied for these models. 
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